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Abstract — We describe a theoretically-motivated algorithm for 
extrapolation of antenna radiation patterns from a small number 
of measurements. This algorithm exploits constraints on the an- 
tenna's underlying design to avoid ambiguities, but is sufficiently 
general to address many different antenna types. A theoretical 
basis for the robustness of this algorithm is developed, and its 
performance is verified in simulation using a number of popular 
antenna designs. 

Index Terms — antenna pattern; antenna extrapolation; design 
space; Whitney embedding 

I. Introduction 

MEASUREMENT of antenna radiation patterns requires 
extensive sampling of transmitted or received signals 
at many angular locations. Scanning a known antenna over 
the antenna under test is the usual approach, but this can be 
time-consuming. We propose a theoretical basis and a practical 
algorithm for using a small number of measurements to 
estimate the unseen and unmeasured portions of the antenna's 
radiation pattern, by exploiting knowledge about the antenna's 
design. The method is robust enough to detect when the 
antenna's design model is incorrect, yet general enough to treat 
many different popular classes of antennas without essential 
changes. 

Although the inverse source problem is unsolvable in full 
generality, it is solvable if the configurations of radiating and 
passive structures are known to be constrained. In particu- 
lar, a solution of the inverse source problem would infer a 
volumetric current distribution from extensive pattern mea- 
surements. Both current distributions and radiation patterns 
are described by infinite dimensional vector spaces: this leads 
to nonuniqueness. The design spaces we propose have the 
distinct advantage that they are finite dimensional descriptions 
of possible antennas, and therefore constrain the inverse source 
problem to a highly overconstrained problem in its most 
ideal incarnation. Since radiation pattern of an antenna varies 
smoothly as its design is adjusted, the well-known Whitney 
embedding theorenjj applies, and indicates that the inverse 
source problem in this restricted setting is solvable. Further, it 
yields bounds on the number of samples one must take from 
the radiation pattern in order to solve for the antenna's design 
in general. 

Inspired by this general theoretical result, we describe a 
general algorithm that exploits the measurements taken of an 
antenna and its associated design pattern in order to estimate 
the pattern at arbitrary directions. In doing so, we carefully 

(1) SRC, inc., 6225 Running Ridge Road, N. Syracuse, NY 13212 USA 
and University of Pennsylvania, 209 S. 33rd Street, Philadelphia, PA, 19104, 
USA e-mail: mrobinson@srcinc.com, robim@math.upenn.edu 

^for a detailed treatment of the necessary differential topological back- 
ground, see Lee (0 



treat the problem of symmetries that arise both in the antenna's 
configuration and in the structure of the measurement sample 
locations. If unchecked, symmetric sampling methods can 
result in ambiguities in the radiation pattern estimates, but 
they are generally easy to avoid. 

II. Historical development 

Recovery of an antenna's excitation from a limited number 
of pattern measurements, and subsequent prediction of addi- 
tional portions of the pattern has been explored by many other 
authors. In its usual incarnation (for instance, as described by 
1 2 1), the antenna is scanned in its near field. After careful 
calibration for mutual coupling effects, near field scanners 
obtain a very detailed model of the antenna's current distri- 
bution. From this current distribution, accurate models of the 
resulting antenna pattern can be predicted. A disadvantage of 
near-field scanning is that sophisticated, expensive equipment 
and careful calibration are required. (It should be noted that 
if the scanner is designed to fit exactly one specific antenna, 
some equipment complexity can be removed O.) In contrast, 
our approach is general enough to handle many different 
antenna types essentially without modification, can use simple 
equipment (a field strength measurement meter suffices), and 
has minimal calibration requirements. Since the measurements 
can be taken in the far field, mutual coupling is essentially 
eliminated. (Our methodology also works in the near field, 
though the advantages are less clear in that case.) 

Another article with a similar processing chain to ours is 
|4|. In this article, the authors use an explicit far-to-near field 
transformation to filter far-field measurements. In particular, 
since they are working exclusively with planar antennas, any 
excitations found to be off of the plane consitute error, and are 
removed. In essence, they manipulate the excitation structure 
of the antenna under the assumption of a known configuration. 
Their algorithm is therefore linear. Our work is an extension 
of this idea, and allows antenna configuration to vary, resulting 
in the need for nonlinear optimization. 

More theoretically, the successful extrapolation of antenna 
patterns relies on the uniqueness of antenna current distribu- 
tions given a subsampling of the pattern. This was an open 
problem for some time, until solved in the negative by Gbur 
in his dissertation |5|. The particular difficulty is that there is 
too much freedom in specifying the current distribution. There 
exist certain distinct current distributions that result in the 
same radiation pattern. These appear to be highly symmetric 
distributions. They are also nongeneric: patterns that are small 
perturbations of the offending ones uniquely specify current 
distributions. 

We argue that the problem becomes solvable if the currents 
are appropriately constrained; we show that this follows using 
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the techniques of differential topology. Because of the gen- 
erality of the argument we follow, it is possible to choose 
many different parametrizations of the current distributions. 
Because it is most congenial to the design of antennas, we 
divide the parameters into two families: those that describe 
the relative magnitude and phase of signals applied to the 
antenna's ports (excitations in what follows), and those that 
describe the relative positions and sizes of reflecting and 
radiating elements (the antenna's configuration). 
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A. Contrasting our approach with previous work 

Unlike the research previously discussed, our approach is 
substantially more general. Since it is borne of a topological 
treatment of the problem, our approach is generically insensi- 
tive to the type of antenna and associated measurement cam- 
paign. Algorithms based on this topological theory of radiation 
patterns are therefore fiexible enough to handle any kind of 
antenna, requiring the authoring of only a pattern computation 
model and a representation of the antenna design space. Our 
algorithms demonstrate good performance on a wide variety 
of interesting antennas, rather than being effective only on a 
particular kind of antenna. Although we frame these examples 
from the point of view of far-field, scalar measurements, 
this is by no means necessary. Indeed, our methods include 
polarization- sensitive measurements and near-field scanning as 
special cases. 

III. General theory 

In this section, we develop the general theory that permits 
extrapolation of antennna patterns. We begin by giving precise 
descriptions of antenna design spaces and antenna patterns, 
then we apply the Whitney embedding theorem to this setting, 
and finally we describe a practical algorithm derived from this 
theoretical grounding. 

Our approach relies on the fact that when enough pattern 
measurements are taken, the radiating structure of the antenna 
is completely constrained. Given this radiating structure, we 
can extrapolate the directivity at any location. Of course, a 
crucial point is that it takes finitely many parameter to specify 
a particular antenna in a class of possible designs. 

A. Design and measurement spaces 

A design space is a product C x £^ of a finite dimensional 
manifold C (possibly with boundary) called the configuration 
and a vector space of complex A^-tuples £ = C^, the 
excitation. We interpret the configuration as defining the pa- 
rameters of an antenna's design: perhaps the element positions, 
locations of reflectors, type of grounding, etc. The excitation 
consists of the relative magnitudes and phases being applied 
to each active element (or feed port) of the antenna. 

Here are some typical examples of antenna designs and their 
associated design spaces: 

1) General phased array with N elements (location of each, 
overlapping elements are permitted): C = R^^, and £ = 

2) Rectangular array (see Figure [T]) with M rows and N 
columns: C = R^+^-2, and 8 = C^^. (Note that the 



Fig. 1. Schematic of configuration space for rectangular phased arrays 
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Fig. 2. Schematic of configuration space for E-plane horn 



number of row spacings is M — 1 and the number of 
column spacings is A^ — 1.) 

3) Horn (see FigureB: C = R^ smdE = C. 

4) Dish (see Figure [3] for a picture of a planar slice of the 
antenna; the full configuration includes the perpendicular 
radius and curvature, as well as orientation of feed). 
Assuming a fixed, known feed with an elliptical dish 
cut from a paraboloid: C = R^ x R^ x R^ >< gQ(^^^ ^^^^ 

£ = C (SO (3) is the group of orthogonal 3x3 matrices 
with unit determinant, and also the group of rotations in 
3 dimensions.) 

The primary focus of this article is an antenna's pattern, 
which we axiomize as a smooth map T from the design space 
P to a space of measurements M, a Banach space of V- 
valued functions on the unit sphere. We assume F is a C- 
vector space, describing the polarization response in the far 
field. We assume that J^ is linear on the £ factor, and may 
be (usually is) nonlinear on C. Generally speaking, V can be 
less informative depending on the kind of pattern under study: 
F = R for a magnitude-only pattern, F = C for a pattern in 
magnitude and phase, V = C^ for a fully-polarized response. 

Of course, J^ is a theoretical construct only; in practice one 
takes measurements of the pattern at finitely many directions. 
Each one of these sample points is an element of {^i}|Li C 



Dish radius 



Configuration 
also includes 
radii in plane 
perpendicular to 
one shown here 



^ Feedpoint location ,' 



Radius of curvature 



Fig. 3. Schematic of configuration space for isotropic-fed parabolic dish 
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S'^. In this case, the full pattern T restricts to the (sampled) 
antenna pattern V :V ^ V^ . 

One should be aware that our model does not explicitly in- 
corporate any physical model of propagation, mutual coupling, 
or excitation. The approach discussed here is topological and 
thereby insensitive to the physics beyond what we assumed 
above. As will be clear in the next section, the critical resource 
for extrapolating patterns is the dimension of the number of 
sample points K versus the dimension of V. 

B. Embedding of the design space 

The obstruction to determining an antenna's design from 
measurements is lack of injectivity of the pattern V. However, 
it is a general fact that V is injective when there are enough 
sample points. 

Theorem 1. Suppose that 2 dim V < K dim V. For generic 
choice ofV^ C^(P;F^), there exists a smooth extrapola- 
tion map To {V)~^ : V{V) -^ M from a submanifold ofV^ 
to M. 

In particular, this map takes a sampled pattern in the image 
of V to the pattern extended to all possible sample point 
locations. This map is well defined and smooth for generic 
choices of V. 

Proof: First, observe that the pattern V : V ^ V^ 
is a topological embedding (hence injective) under a generic 
perturbation in C^(P; V^) when 2dim7:> < K d\mV . This 
follows immediately from the Whitney embedding theorem as 
given in |6| for manifolds with boundary, since 7^ is a smooth 
map. One should notice that although we assumed smoothness 
of V (dependence of the pattern on the configuration), in fact 
we have obtained smoothness of the inverse. 

It then follows immediately that J^ is generically an embed- 
ding, since M is an infinite dimensional vector space. ■ 

The reader is cautioned that genericity of V is not a 
completely innocuous condition. Precisely, it indicates that in 
the vast majority of situations, patterns will uniquely determine 
antenna configuration and excitation. However, this implies 
that in the situations where V is injective, there is some 
measure of symmetry-breaking. 

A simple example of a choice of V for which the conclusion 
of Theorem [T] does not hold arises from sampling azimuth 
patterns only. Specifically, let us consider the case of a 
phased array composed of N isotropic radiators located at 
{xm}m^i with excitations {am}m^i respectively. Consider 
the pattern map arising from taking sample points at 5^ = 
(i^cos kO^Rsin kO^ 0) (where R^ are fixed), namely 

N 



V{xi, ...,XAr,ai, ...,aN)k = ^ <^r) 



There is an inherent ambiguity whether a given Xm is above 
and below the x?/-plane. On the other hand, many small 
vertical perturbations of the sample points will induce a 
small perturbation of the resulting pattern maps (since exp 
is smooth), yet will break the symmetry and thereby ensure 
the existence of a extrapolation map. (Proving injectivity of 



V in this case is an elementary, if tedious, exercise for the 
reader.) 

Because symmetric (or nearly symmetric) sampling is com- 
mon in practice, the existence of a extrapolation map could 
be an inherently delicate phenomenon. In order to validate 
its robustness in practice, extensive simulation is discussed in 
Section |IV| with realistic antenna configurations and various 
sampling patterns. 

C Algorithm for extrapolating the pattern 

In this section, we describe an algorithrrj^for computing and 
approximation to the inversion V~^ : V{V) C V^ -^ V. In 
this way, we determine the configuration and excitation of the 
antenna from its pattern. Since V may have a large dimension, 
and V may have a complicated form, it is desirable to devise a 
method that is flexible, yet exploits structure in 'P. To this end, 
we use the product decomposition V = C x £, and the fact 
that V depends linearly on the £ factor. This has two major 
benefits: (1) constrained least-squares can be used to robustly 
recover the element of 8 given a known element of C, and (2) 
the nonlinearity of V is confined to a much lower dimensional 
factor, namely C. We can therefore use a gradient-like search 
algorithm to traverse C in a robust fashion. 

Our algorithm takes as input 

• A parametric representation of the antenna design space 
V, specifically written as the product C x £, 

• The locations of each sample point, 

• Signal measurements at each sample point, 

• A sample point for which an estimated measurement is 
desired, and 

• An implementation of the pattern predictor T. 

As output, it generates an estimate of the measurement at 
desired sample point. 

The algorithm proceeds recursively, starting with (xq, ao) G 
C X S and an associated set of measurements p e M, and 
producing a sequence (xn^an) which we compute until some 
convergence criterion is met. Each iteration is computed from 
the iteration immediately before it, in the following way. 

1) A dense rectangular grid with a fixed number of el- 
ements of C is computed, called {xi^}j, centered on 
Xn- This grid's orientation depends on the particular 
parametrization of C that was chosen for the antenna 
under test. The spacing between adjacent x{^ elements 
is chosen to be a decreasing function of n. 

2) For each element of {x:^}, the best (in the least squares 
sense, constrained to unit magnitude) excitation a^^ asso- 
ciated to the linear problem T{x{^, a^) = pis computed. 

3) The next iterate (xn+i,an+i) is defined to be 
argmin^^j^^^j^^ || J'(4, a{) - p||. 

IV. Simulation 

We have validated performance of our algorithm using 
appropriately parametrized design spaces V, as described in 
the previous section. These results confirm both the theoretical 
and practical value of our method, by showing that the pattern 

^A patent is pending on the algorithm 
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Fig. 4. Spatial layout of sampling patterns used in our study, clockwise 
from top left: compete azimuth and elevation samples, blocks of contiguous 
azimuth samples, random samples. Axes are in meters. Not shown: azimuth 
samples only 



Fig. 5. A typical evolution of residuals over solver iterations for a 3x3 
rectangular phased array. Varying amounts of measurement error was added 
to the sampled patterns to test robustness of convergence. 



map V is injective in these cases, and can be used to robustly 
extrapolate to an approximation to T. 



A. Sampling patterns 

In order to determine sensitivity of our algorithm to the 
choice of sample points, we tried several kinds of sample point 
layouts. Specifically, we examined (see Figure |4]) 

1) azimuth pattern samples only, 

2) blocks of contiguous portions of azimuth samples, at 
random small elevation angles, 

3) complete azimuth and elevation patterns together, 

4) and random sampling. 

It's worth noting that although azimuth and elevation patterns 
are considered standard procedure due to the mechanical 
simplicity of collecting them, they can be suboptimal for our 
techniques. Specifically, if the azimuth- or elevation-plane is 
a plane of symmetry for the antenna, ambiguities in T may 
arise. Therefore, to balance between measurement complexity 
and the need to break sampling symmetry, we also considered 
contiguous blocks of azimuth samples at random elevation 
angles. We compared the performance of our approach using 
a large number of typical antennas to uncover the relationship 
between sampling patterns and the resulting ambiguities. 

B. Performance metrics 

Our experimental procedure aims to validate the perfor- 
mance of the extrapolation algorithm by comparing its output 
to simulated truth. We did this by generating a densely sampled 
pattern at many azimuth and elevation angles from a known 
configuration. From this densely sampled pattern, the subset 
of samples to be used as input to the algorithm was selected. 
The algorithm was then applied to estimate the pattern at each 
of the original sample points. Although directly computing 
the differences between the estimated and truth values gives a 
measure of performance, the total error is not very useful in 
practice, since the truth pattern is unavailable. 



A more predictive performance measure exploits the fact 
that the estimated pattern differs from the suppled (under- 
sampled) pattern. This pattern residual error (between output 
of the algorithm and the input pattern) should ideally be 
zero, but due to discretization and numerical errors does not 
vanish. By comparing the residual error and the total error, it 
becomes possible to predict the performance of the algorithm 
knowing only the residual error. Indeed, in our example cases, 
reliable estimates can be obtained. For instance, consider 



Figures |7| |9] and [T0| As the residual (the difference between 
supplied pattern and the predicted pattern at the sample points) 
decreases (moving leftward in Figures [7] |9] and[TO]), it is clear 
that the total error in the pattern generally decreases. 

It should be noted that the antenna pattern measurement 
space V^ can be a vector space over C (magnitude and 
phase) or R (magnitude). In the latter case, it is theoretically 
immaterial whether the magnitude is given in absolute units 
or in decibels. 

Additionally, one can examine the residual errors and their 
evolution over the course of the iterations of the algorithm. In 
Figure |5j several different amounts of additive Gaussian white 
noise was added to the pattern measurements, and the resulting 
degradation of residual error was recorded. Since the residual 
error decreases over the course of the algorithm's execution 
without much dependence on measurement error, we deem the 
algorithm to be robust to these errors. 

Residual error evolution can also be used to identify in- 
consistences in the choice of design space and thereby make 
refinements in the model. This permits substantially greater 
power than is indicated by the theory. For instance, if a 
phased array is correctly postulated for the unknown antenna's 
configuration, but the wrong number of elements is specified, 
it will be difficult to reduce the residual to zero. Supplying the 
correct number of elements permits the algorithm to converge. 

V. Results 

In this section, we discuss a number of special cases of 
Theorem [T] We indicate what the configuration space is for 
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the particular antenna in question, as well as parametrizations 
we have found convenient for use in our algorithm. 

To keep the exposition simpler, we treat idealized models of 
antennas, ignoring polarization, ground effects, support struc- 
tures, and dielectrics. The reader may supply more realistic 
models as appropriate. 

A. Rectangular phased arrays 

As noted in Section |III-A| rectangular phased arrays have 
a design space given by C = R^+^-^ and S = C^^, 
where there are M rows and N columns of elements. The 
components of C specify the element spacing, as shown in 
Figure [T] In this case, the full pattern map is straightforward 
to construct, namely 

M N 



J='{x, y, a){(t), 0) = ^^ ajkC 



ioj{yj sin 9-\-Xk sin 0)/c 



3 = 1 k=l 

where Xk represents the column locations (xi = is fixed), i/j 
represents the row locations (yi = is fixed), and the antenna 
elements lie in the x?/-plane. 

It is known that the full pattern of a given rectangular 
phased array is uniquely determined by its excitations. This 
follows from the fact that the near-to-far transform (taking 
£ -^ L^(S'^,C)) is a discrete Fourier transform under an 
appropriate change of variables, and the Fourier transform is 
invertible. (See 1 7 1 for a standard treatment.) Our Theorem [T] 
extends this standard result to treat subsampled patterns, and 
our algorithm can treat this situation. Generically, one would 
need to take K measurements, where K > 2 -\- 2MN. 

To validate performance in simulation, we tested our al- 
gorithm on small arrays with 3 rows and 3 columns of ele- 
ments, spaced randomly up to one-half wavelength. Random 
complex excitations were applied to each element. A typical 
antenna pattern, its principal plane patterns, and extrapolation 
results are shown in Figure [6] Typical performance of our 
algorithm (using 10 iterations on 100 sample antennas) on 
these rectangular arrays is shown in Figure [7] It is clear that 
the algorithm works well: if the residual error is low, the total 
error will usually be low also. However, there is substantial 
dependency on the sampling pattern. Principal plane patterns 
do not exhibit a definite reduction in total error, even with 
essentially no residual errors. This is due to the vertical and 
horizontal symmetry of the array, which leads to ambiguities 
in the resulting pattern when sampled in the azimuthal and 
elevation planes. Clearly, a lower bound on total error is 
achieved with random sampling. 

It therefore is important to be aware of these symmetries 
when selecting collections of sample points. In particular, 
using strictly azimuth or elevation patterns (or both) can result 
in mis-estimation of non-cardinal sidelobes as shown in Figure 
[8] These sidelobe mis-estimations are the cause of the large 
total errors in Figure [7] We therefore recommend using a 
different sampling pattern for rectangular arrays to avoid these 
issues. 

Rectangular arrays present an additional challenge: the 
number of elements may be unknown. However, in this case, 
the residual errors cannot be decreased beyond a certain 
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Fig. 6. Typical rectangular array pattern. Top frames show azimuth and 
elevation patterns respectively; dotted lines represent the closest match found 
by our algorithm. Bottom left frame shows truth pattern, bottom right is the 
prediction from our algorithm. 
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Fig. 7. Extrapolation performance for rectangular arrays. Accumulation of 
points near vertical axis indicates the presence of persistent ambiguities with 
certain sampling patterns. 
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Fig. 8. Incorrectly predicted off-axis sidelobe levels for a rectangular array, 
due to symmetric sampling. Top frames show azimuth and elevation patterns 
respectively; dotted lines represent the closest match found by our algorithm. 
Bottom left frame shows truth pattern, bottom right is the prediction from our 
algorithm. 
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TABLE I 

Value of smallest residual given varying numbers of rows and 
columns of elements (truth antenna has 3 columns and 4 

ROWS) 



Total versus residual error (Horn antennas) 





1 column 


2 columns 


3 columns 


4 columns 


5 columns 


1 row 


730 


469 


283 


154 


145 


2 rows 


418 


261 


47 


47 


47 


3 rows 


397 


225 


2 


0.5 


0.5 


4 rows 


390 


224 


0.5 


0.4 


0.3 


5 rows 


360 


216 


0.4 


0.3 


0.2 



amount. This provides a discrimination mechanism: if the 
residual errors do not decrease given a number of iterations 
and random initial conditions, a different number of elements 
in the configuration should be tried. To test this methodology, 
we simulated an antenna with 3 columns and 4 rows of 
elements. We then computed the minimal residuals with our 
algorithm for all configurations with between 1 and 5 rows and 
columns. These residuals are shown in Table ID It should be 
clear that using more rows in the proposed configuration than 
are actually present will result in low residuals. Therefore, 
the correct configuration is the one whose residuals are the 
smallest and has the smallest number of rows and columns. 
In our table, this occurs at precisely the correct number of 
elments (3 columns, 4 rows). 

B. Horn antennas 

There are several different kinds of horn antennas that are 
in typical usage. For concreteness, consider an E-plane horn. 
This kind of antenna has three design parameters (see Figure 
[2]): the width between parallel plates, the height of the throat of 
the horn, and the height of the mouth of the horn. Therefore, 
C = M^. In contrast to a phased array, the excitation of a horn 
antenna presents only one complex degree of freedom. 

From these parameters, it is straightforward to compute the 
pattern J^ of the horn. For our simulation, we used the closed- 
form solutions given by (13-1 lb) or (13-1 Ic) in lO. For this 
antenna. Theorem [T] indicates that one needs to take at least 10 
measurements to obtain injective pattern maps V. In practice, 
errors will be present so more measurements are generally 
desirable. 

The horns in our simulations are oriented so that the E- 
field is in the elevation plane. Typical horns with mouth 
and side length dimensions up to 6 wavelengths were tested. 
Overall performance of the extrapolation of E-plane horn 
antennas is shown in Figure [9j where 100 antennas were 
tested and the solver used 10 iterations at most. Unlike the 
case of rectangular phased arrays, there do not appear to be 
any prominent symmetries that degrade performance when 
principal planes are used. However, using contiguous blocks of 
azimuth samples (with random small elevation angles) appears 
result in some ambiguity. Again, random sampling (which 
has the least amount of symmetry) acheives a practical lower 
bound on total error. 

C Dish antennas 

Dish antennas can be complicated by the fact that the 
illumination of the dish itself is intimately connected to the 
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Fig. 9. Extrapolation performance for E-plane horns. Accumulation of points 
near vertical axis predicted with azimuth blocks indicates the presence of an 
ambiguity with this sampling pattern. 



feed structure. As explained in Section III-A C = M x R^ x 
R^ X 6*0 (3) (or similar) in this case, as shown in Figure [3] 
Unlike the case of phased arrays or horn antennas, no closed 
form solutions typically exist. Therefore, the computation 
of J^ in our algorithm will usually involve a computational 
electromagnetics engine. In any case, one needs to make at 
least 20 measurements to obtain an injective pattern map. 

To make this model concrete enough for a computer simu- 
lation, we represented the feed as an isotropic radiator. In this 
case, the S0{?>) factor plays no role, and may be removed 
from C. Indeed, the radiation pattern of the resulting isotropic- 
fed dish can be found by ray-tracing from the feed to the 
receiver, reflecting off each of a dense set of points on the 
reflector. Specfically, the reflector is described by the equation 
of y = ax^ + hz^. 

Our simulation tested dish antennas whose reflector radii 
were up to 6 wavelengths and feed was located within 6 
wavelengths from the vertex of the reflector. The parameters 
a and h were chosen randomly between and 6. The reflector 
was simulated with a grid with 20 axial points and 10 radial 
points. The typical results for 100 dish antennas (10 iterations 



used in the solver) are shown in Figure 10 Unlike both horn 
and rectangular arrays, the performance of dishes is essentially 
the same regardless of sampling pattern. In particular, random 
sampling does not outperform other sampling patterns. There 
are no obvious ambiguities, apparently as the model of the 
antenna has no symmetries that always are shared the sam- 
pling pattern. The model can be symmetric, but pairs typical 
antennas will not share axes of symmetry. 

VI. Conclusion 

This article presents a theoretical and algorithmic frame- 
work for extrapolating antenna radiation patterns from sub- 
sampled patterns. Notionally, we aim to dramatically improve 
the speed of measurement of full patterns taken from a far- 
field range, allowing off-cardinal axis sidelobe structure to be 
estimated indirectly. Our approach is based on the Whitney 
embedding theorem of differential topology, and is therefore 
sufficiently general to treat many different kinds of antennas 
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Fig. 10. Extrapolation performance for isotropic fed parabolic dishes. No 
persistent ambiguities are present in the pattern predictions. 



with no change in algorithms. We vaUdated the performance 
of our algorithms in simulation for three different kinds 
of antennas: rectangular phased arrays, E-plane horns, and 
parabolic dishes. 
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